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The  method  of  multiple  scales  is  used  to  determine  a  first-order 
uniform  expansion  for  the  effect  of  counter-rotating  steady  streamwise 
vortices  in  growing  boundary  layers  on  Tollmien-Schlichting  waves.  The 
results  show  that  such  vortices  have  a  strong  tendency  to  amplify 
three-dimensional  Tollmien-Schlichting  waves  having  a  spanwise  wave¬ 
length  that  is  twice  the  wavelength  of  the  vortices.  An  analytical 
expression  is  derived  for  the  growth  rates  of  these  waves.  These  growth 
rates  increase  linearly  with  increasing  amplitudes  of  the  vortices. 

I.  INTRODUCTION 

We  consider  the  effect  of  counter-rotating  steady  streamwise 
vortices  on  the  instability  of  growing  boundary  layers.  We  describe 
a  parametric  instability  mechanism  by  which  such  vorticies  amplify 
selected  three-dimensional  Tollmien-Schlichting  waves.  To  first  order, 
the  selected  waves  have  a  spanwise  wavelength  that  is  twice  that  of  the 
vortices. 

Weak  and  moderately  strong  steady  streamwise  vortices  arise  abun¬ 
dantly  in  boundary  layers  from  many  causes.  In  a  series  of  wind-tunnel 
tests  over  flat  plates  with  zero-pressure  gradient,  Klebanoff  &  Tidstrom 
(1959)  observed  steady  quasi-periodic  variations  in  the  spanwise  direc¬ 
tion  (streamwise  vortices)  evidently  evoked  by  freestream  conditions. 
Similar  vortices  were  observed  in  a  National  Physical  Laboratory  tunnel 


specifically  designed  for  the  study  of  two-dimensional  boundary  layers. 
Bradshaw  (1965)  found  that  these  variations  may  appear  downstream  of 
slightly  nonuniform  settling-chamber  damping  screens,  depending  on  their 
solidity.  Using  the  method  of  matched  asymptotic  expansions.  Crow  (1966) 
inferred  the  effect  of  a  small,  periodic  incident  transverse  flow  on  the 
mean  boundary  layer  over  a  flat  plate. 

Gortler  (1941)  found  that  a  boundary  layer  over  a  concave  surface 
is  strongly  unstable.  The  instability  is  manifested  by  the  presence  of 
counter-rotating  vortices  having  their  axis  in  the  streamwise  direction. 

Using  the  tellurium  method,  Wortman  (1969)  gave  detailed  flow  visualization 
of  these  vortices  on  slightly  curved  walls.  Bippes  ( i 973 )  conducted 
experiments  on  walls  with  radii  of  curvature  of  0.5  and  1m  so  that  the 
generated  Gortler  vortices  are  fairly  strong.  He  made  the  flow  visible 
by  using  the  hydrogen-bubble  technique  and  photographed  it  with  a  photo- 
grogrammatic  sterocamera.  He  analyzed  photogranmetrically  the  photographs 
and  obtained  fairly  accurate  quantitative  information  of  the  flow  field. 

Unlike  the  case  of  preexisting  streamwise  vo*tices,  Gortler  vortices  generated 
by  a  concave  surface  are  amplified  with  streamwise  distance.  Their 
amplification  is  exponential  (Smith,  1955)  when  they  are  weak  and  it 
appears  to  be  linear  when  they  are  strong  (Bippes,  1978).  Floryan  and 
Saric  (1979)  gave  a  comprehensive  review  of  the  different  analyses  of 
these  vortices  and,  using  the  model  of  Smith  (1955),  presented  fairly 
accurate  numerical  results  describing  these  vortices. 

In  their  experiments  on  the  nature  of  boundary- layer  stability, 

Klebanoff,  Tidstrom  &  Sargent  (1962)  established  that  streamwise  vortices 
are  associated  with  nonlinear  three-dimensional  wave  motions.  Benney  & 

Lin  (1966)  modeled  the  generation  of  these  vortices  by  the  nonlinear 
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interaction  of  a  two-dimensional  wave  with  a  three-dimensional  wave 
superimposed  on  the  laminar  profile.  Using  this  model,  Antar  &  Collins 
(1975)  calculated  these  vortices  in  a  boundary  layer  on  a  flat  plate 
for  different  amplitudes  of  the  waves. 

A  number  of  experimental  studies  investigated  the  influence  of 
steady  streamwise  vortices  on  the  transition  from  laminar  to  turbulent 
flow.  Aihara  (1962),  Tani  &  Sakagami  (1964),  and  Tani  &  Aihara  (1969) 
studied  the  influence  of  steady  streamwise  vortices  on  two-dimensional 
Tollmien-Schlichtii.g  waves  generated  by  a  vibrating  ribbon.  They  con¬ 
sidered  the  case  in  which  the  vortices  were  generated  naturally  on  a 
concave  surface  (GSrtler  vortices)  as  well  as  the  case  in  which  the 
vortices  were  generated  artifically  by  a  row  of  wings  on  a  cancave 
surface.  They  measured  the  distributions  of  the  mean  velocity  and  wave 
intensity  across  the  boundary  layer  for  three  spanwise  positions  at  a 
number  of  streamwise  stations.  They  concluded  that  the  Gortler  vortices 
indirectly  affect  the  transition  by  inducing  a  spanwise  variation  in 
boundary-layer  thickness,  at  least  when  the  radii  of  curvature  are  not 
extremely  small.  However,  they  did  not  present  any  measurements  of 
the  growth  rates  of  the  Tollmien-Schl ichting  waves.  Wortman  (1969) 
investigated  the  development  of  natural  transition  downstream  of  Gortler 
vortices.  Using  the  tellurium  method,  he  determined  the  direction  and 
relative  magnitude  of  the  unsteady  velocities  from  the  streaklines  by 
confining  his  observations  to  the  vicinity  of  the  starting  point  of  the 
streaklines.  He  observed  a  steady  second-order  instability  that  destroys 
the  symmetry  of  the  Gb'rtler  vortices.  He  suggested  that  this  instability 
is  caused  by  secondary  vortices  having  spanwise  wavelengths  that  are 


twice  those  of  the  Gfirtler  vortices.  Then  he  observed  a  third-order 
instability,  consisting  of  regular  three-dimensional  oscillations. 

The  above  shows  that  there  are  many  theoretical  and  experimental 
studies  relating  to  the  generation  of  streamwise  vortices  and  a  number 
of  experimental  studies  relating  to  their  effect  on  transition,  but  to 
the  author's  knowledge,  no  theory  yet  exists  on  how  these  vortices 
affect  the  development  of  Tollmien-Schlichting  waves.  The  purpose  of 
the  present  paper  is  to  present  a  parameteric  instability  mechanism  by 
which  the  streamwise  vortices  increase  the  growth  of  selected  Tollmien- 
Schlichting  waves  in  growing  boundary  layers.  To  first  order  the 
selected  waves  have  a  spanwise  wavelength  that  is  twice  that  of  the 
vortices,  while  to  second  order  the  selected  waves  have  spanwise  wave¬ 
lengths  that  are  equal  and  twice  that  of  the  vortices.  To  minimize  the 
algebra,  we  consider  flows  with  growing  steady  counter-rotating  vortices 
over  flat  plates  with  or  without  pressure  gradients. 


II.  Problem  Formulation 


We  consider  the  stability  of  a  basic  flow  that  consists  of  the 
superposition  of  the  Blasius  or  Falkner-Skan  flow  and  a  flow  corres¬ 
ponding  to  growing  steady  quasi-periodic  counter-rotating  streamwise 
vortices.  Thus  we  consider  the  stability  of  the  flow  described  by 


U  =  U0(xi,y)  +  Cylii  (xi  ,Zi  ,y)cos2Bz  +  ....  (1) 
V  =  eVo(x2 ,y)  +  evVi(xi ,Zi ,y)cos2Bz  +  ....  (2) 
W  =  evWi(xj,Zi,y}sin2Bz  +  ...,  (3) 
P  =  P0(xi)  +  evPi(xi,Zi,y)cos2$z  +  ....  (4) 


where  the  subscript  o  refers  to  the  Blasius  of  Falkner-Skan  flow,  the 
subscript  1  refers  to  the  flow  corresponding  the  the  streamwise  vortices, 
xi  =  ex,  Zi  =  ez,  e  is  a  small  dimensonless  quantity  that  accounts  for 
the  quasi-periodicity  of  the  vortices  as  well  as  for  the  streamwise  growth 
of  the  vortices  and  the  boundary  layer,  S  is  a  real  dimensionless  spanwise 
wavenumber,  and  ey  is  a  small  dimensionless  quantity  that  indicates  the 
strength  of  the  streamwise  vortices.  In  Eqs.  (1)  -  (4),  velocities  and 
lengths  are  made  dimensionless  by  using  a  reference  velocity  Ur  and  a 
reference  length  6r<  We  superpose  the  small  unsteady  perturbation 
quantities  £yu(x,y,z,t) ,  £yv(x,y,z,t) ,  eTw(x,y,z,t) ,  and  Cyp(x,y,z,t) 
on  those  given  in  Eqs.  (1)  -  (4)  so  that  the  total  flow  quantities 
become  U  +  CyU,  V  +  CyV,  W  +  CyW,  and  P  +  Cyp.  Here,  Cy  is  a  small 
dimensionless  quantity  that  is  the  order  of  the  amplitude  of  the  Tollmien- 
Schlichting  waves.  In  this  paper,  Cy  is  assumed  to  be  much  smaller  than 
ev  and  e  so  that  terms  the  order  of  c|  can  be  neglected  compared  with  CyCv 
and  eye.  Substituting  these  total  flow  quantities  into  the  dimension- 


less  Navier-Stokes  equations,  subtracting  the  basic-flow  quantities,  and 
keeping  linear  terms  in  the  perturbation  quantities,  we  obtain 


3u  3v  8w  _  n 
3x  3y  3z  '  ’ 


If  *  «  If  *  »  37  +  »  0  +  *  |  *  “  ST  *  -  If  ■  -  If  *  “ • « 5 1 

If  +  u  If  +  u  H  +  v  0  +  v  I?  *  “  If  *  w  ll  ■  -  -  ^  v  • ' 7 > 


3y  R 


|wtu|S*u|Stv|S*vf*w||t>.|S=-||*l  V*«,<8) 

where  t  is  made  dimensionless  by  using  5r/Ur,  R  =  Ur6r/v  is  the  Reynolds 
number,  and  v  is  the  kinematic  viscosity  of  the  fluid. 

Substituting  Eqs.  ( 1 )- (4)  into  Eqs.  (6)-(8),  we  obtain 

If  *  «•  If  * » &  ♦  If  -  i  72“  *  %  If  ♦  »>  0 ♦  It  ’']cos2ez 

*7  ^ 

♦  [Mi  H  -  2eU,w]s1n2fSz  -  e[|^  u  +  V,  |^]  +  0(eev)  (9) 

If +  u*  0 4  0  - 1 7!v  •  -  Ev|i°>  If  *  "■  0 +  It  v]cos28z 

-»  - 


+  [Wi  |y  -  2BViw]sin28z 


On.  +  II n  —  +  -  1  vzw  =  -  e 

3t  Uo  3x  3z  R  w  V 


+  [|£-  V  +  W,  I^]sin26z  \  -  e  V„  f=-  +  0(eew) 


zj  -  £tv»  |y  +  |y- 


v]  +  0(e2,eev) 


+  Vi  §=•  +  2f5Wiw]cos28z 


Equations  (5)  and  (9)-(ll)  need  to  be  supplemented  by  initial  and  bound¬ 
ary  conditions.  The  initial  conditions  are  specified  later,  while  the 
boundary  conditions  for  an  impermeable  flat  surface  are 

u  =  v  =  w  =  0aty  =  0,  (12) 

u,  v,  w  -*■  0  at  y  -*•  °°.  (13) 
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III.  Solution 


We  use  the  method  of  multiple  scales  (e.g.,  Nayfeh,  1973)  to 
determine  a  first-order  uniform  expansion  for  Eqs.  (5)  and  (9)-(13).  To 
accomplish  this,  we  let  ey  =  0(e)  and  write  ey  =  X£.  where  x  =  0(1).  If 
c  <<  e  the  effect  of  the  growth  of  the  boundary  layer  is  small  com¬ 
pared  with  the  effect  of  the  vortices.  If  ev  «  e.  The  effect  of  the 
vortices  is  small  compared  with  that  due  to  the  growth  of  the  boundary 
layer,  and  the  solution  accounts  for  the  nonparallel  effects  only.  Thus 
the  above  ordering  yields  an  expansion  that  accounts  for  the  effects  of 
the  streamwise  vortices  and  the  growth  of  the  boundary  layer,  and  it 
includes  the  cases  e  «  ey  and  cy  «  c  as  special  cases. 

We  seek  a  uniform  expansion  for  Eqs.  (5)  and  ( 9) - ( 1 3 )  in  the  form 

u  =  {  enu_(xo »Xi ,y,z0 ,zi ,t0 ,ti )  +  0(e2), 

n=0  " 

v  =  J  env  (xo ,Xj ,y,Zo ,Zi ,to ,ti )  +  0(e  ), 
n=0 

w  =  J  e^w  (xo >Xi ,y,Zo ,Zi ,to ,ti )  +  0(e  ), 
n=0 

p=  l  enpn(xo,Xi,y,z0,z1,t0,ti)  +  0(c2), 

n=0 

where 

y  =  rnx  z  =  Cnz,  t  =  ent  . 

xn  E  x’  Zn  n 

Substituting  Eqs.  (14)-(18)  into  Eqs.  (5)  and  (9)-(13)  and  equating 
coefficients  of  like  powers  e,  we  obtain 


(14) 

(15) 

(16) 

(17) 

(18) 


7 


Order  e° 


^(ut.Vo.p,)  Vju,  ■  0, 

-Tstvo.Po)  •  lit  +  U,  g-  t  1^-  -  1  7i»,  =  0, 

=  o, 

u0  =  v0  =  w0  =  0  at  y  =  0, 
u0,v0,w0  -*-0  as  y 


09) 

(20) 

(21) 

(22) 

(23) 

(24) 


Order  e 


/l(u, =  -  lif  -  Hf  . 

-  x^u.  &  *  ».  If  *  f1-  V.1COS26Z.  ♦  [W,  |ax 

-  2BUiW(]sin2BzJ  -  u0  -  V,  |jl  , 

/,(v„p,)  -  -  &  -  u.  &  *  |  [077  +  life;!  -  xjiu. 

+  Vl  ly^  +  fy1-  yo3cos2Bzo  +  [H,  ||j-  -  2gV,w0]sin26zoj 

-Vo|^-^v0, 

*<"■  ’P* >  =  *  f?7  -  u0  +  I  fe  +  fe] 


(25) 


(26) 


3v0 

3x0 


(27) 


-  xtfu,  f*J-  +  V1  +  28Wiw0]cos2Bz0  +  [f^  v0 


9w 


-3W, 


+  “■  0si"i»2]  -  V*  If  • 


(28) 
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u i  =  v i  =  w i  =  0  at  y  =  0, 
Ui.Vi.W!  0  as  y 


(29) 

(30) 


where 


Vo  = 


?  ^2  d? 


72  _  3'  .3 


3Xo  3y  3Zo 


The  initial  conditions  are  taken  such  that  the  solution  of  the 
zeroth-order  problem,  Eqs.  ( 1 9 ) - ( 24 ) »  consists  of  two  wave  packets 
centered  around  the  frequency  w,  the  streamwise  wavenumber  a,  and  the 
spanwise  wavenumbers  8i  and  -  gi ;  that  is, 

u0  =  Aj (x, ,z, ,tj )c, ) (y,x, )exp(ie, )  +  A2(xI,z1,t1)521(y.x1)exp(ie2), 


Vo  =■  A1c12(y.x1)exp(i6i)  +  A2c22  (y,x3  )exp(ie2 ) , 


w0  =  A1r,13(y,x1)exp(i6i)  +  A2^2  3  (y,x3  )exp(i02) , 


Po  =  A1Cm(y,x1)exp(i61)  +  A2c2 „ (y ,x3  )exp(ie2 ) , 


where 

Gn  =y 


0„  =  | ot(xj  )dx  -  Bnz0  -  cuto,  82  =  -  8i  » 


(31) 

(32) 

(33) 

(34) 

(35) 


and  the  functions  Ai  and  A2  are  undetermined  at  this  level  of  approxima¬ 
tion;  they  are  determined  by  imposing  the  solvability  conditions  at  the 
next  level  of  approximation.  Substituting  Eqs.  (31 ) - ( 35 )  into  Eqs. 

( 1 9)- ( 24 )  yields  the  following  eigenvalue  problems: 
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+  °Si2  -  i6n^n3  =  0  * 


i(U0a  -  wkni  +  e^DUo  +  iar,n4  -  ^  (D2  -  a2  -  B?)cni  =  0,  (37) 

i(U0a  -  w)cn2  +  Dcni<  -  ^  (D2  -  a2  -  8?)^  =  0,  (38) 

Huoa  -  „Knj  -  feni„i<  -  J  (d2  -  -  e?kn)  ■  o,  (39) 


=  ?n  =  Cn  =  0 

n  i  n  2  n  3 


WS.3  - 0 


at  y  =  0, 


where  Dc  =  9r,/3y.  For  a  given  ui,  8i.  and  k,  one  can  solve  Eqs.  (36)- 

(41)  numerically  to  determine  the  complex  eigenvalue  a  and  the  eigenr 

functions  r  . 

nm 

Substituting  Eqs.  (31 )-(35)  into  Eqs.  (25 ) - ( 28 )  yields 

2  3A  3A 

U  1(U|.V,.W.)  ■  -  l  {—■  c  +  3^  ?„,>«?(«„> 

n=  I 

■  i  AnSTexp(i6n)-  (42) 


r^2(ui,Vi,p1)  -  -  ^  C3ti  cni  +  (uoCni  +  ^n4  R  cnj 


)  _JL 
1  9xi 


+  2iV„,  s>>=<i9n>  -  j,  t(u.  -  ^  -  K,  at;1 

♦  5^  +  'n,  +  V"  5r1]Vx|>(1en)  -  1  xOH.c,;,, 

+  V i  Ci2  -  WiBiCn  +  2i eUlCi  3]A1exp[ie2 

+  2 i  ( S  -  3i  )z o ]  -  \  x[iuiaC2i  +  Vi  l^1-  +  ly1-  C22 
■  Wi3i^2i  ■  2  i  Si)  1 C2  3  ]A2exp[  i  0 1  -  2  i  ( 0  -  Bi)zo]  +  NST,  (43) 
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X 


c5<f  3  (Vi  ,Pl  )  =  -  4n?  +  (U 


,.  3A„  2 1 3  3A 

gigtw  n  _ n  _ n-i 

0  '  R  '  "n2  3xi  R  Si2  3ziJ 


exP(ien)  '  ^  [(Uo  -  ?f)  377 

+  |y  (VoCn?)]exp(iOn)  -  \  xtiaUjCia  +  |y(VxCi2)  -  BjWiC 
+  2i6VirJi3]Aiexp[i02  +  2i(B  -  Bi)z0]  -  j  x[iaUiC22 
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+  yy  (Vi£22)  -  S i Wi ^ 2 2  -  2i BV i ^ 2 3 ]A2exp[ i 0 1  -  2i(B  -  3i)zoj 
+  NST,  (44) 


o^4(wi  ,pi ) 


?  3A 

^  <n,  +  (U“ 


0.  9A 

—  il  +  (r 

R  'c'n3  3xi  i;n, 


2  i  B 

+  C  ) 


R 


n3 ' 


9An-,  /ic  \  v  r/ii  2ia,  ^n3  i  da 
gZi]exp(i6n/  -  [(Uo  -  R  )  3xi  R  r"n3  dxx 


+  Vo  3yJ’^Anexp^Gn^  ’  2~x[iUia^13  +  *  0i)WiCi3 

+  Vi  -  i  -1  Ci2]Aiexp[iG2  +  2i(3  -  Bi)z0] 

-  xC i U 1CXC2 3  +  (23  -  Bi)WiC23  +  Vi  g-y2-3 

+  i  r»2?]A2exp[i0i  -  2i(B  -  Bi)zo]»  (45) 


where  NST  stands  for  terms  that  are  proportional  to  exp[±i(3  +  3i)zo]> 
which  do  not  produce  secular  terms  in  ui,  Vi,  Wi,  and  pi. 

Since  the  homogeneous  parts  of  Eqs.  ( 42 ) - ( 45 ) ,  (29),  and  (30)  are 
the  same  as  Eqs.  ( 1 9) - (24 )  and  since  the  latter  have  a  nontrivial  solu¬ 
tion,  the  inhomogeneous  equations  (42)-(45),  (29),  and  (30)  have  a 
solution  if,  and  only  if,  the  inhomogeneous  parts  are  orthogonal  to 
every  solution  of  the  adjoint  homogeneous  problem.  These  solvability 
conditions  depend  on  whether  S  ~  3i  or  not.  If  3  is  away  from  3i> 
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the  solvability  conditions  yield  two  uncoupled  equations  describing  the 
effect  of  nonparallel ism  on  Aj  and  A2.  If  8  ~  81,  we  introduce  a 
detuning  parameter  o  defined  by 

0  =  Si  +  eo  (46) 


where  a  =  0(1)  and  express  (0  -  20i)zo  as  ozi.  Then,  imposing  the 
solvability  condition  that  the  inhomogeneities  be  orthogonal  to  every 
solution  of  the  adjoint  homogeneous  problem,  we  obtain 

911  fl7  +  912  1x7  +  913  Izf  =  hlAl  +  xhi2A2exp(-iazi)  (47) 

921  Ilf  +  922  Ixf  +  923  fz7  =  hzA2  +  Xh2iAiexp(iaZi)  (48) 


where  the  g's  and  h's  are  given  in  Appendix  A  together  with  the  adjoint 
problems.  Differentiating  Eqs.  (19)  -  (24),  respectively,  with  respect 
to  a  and  Bn  and  imposing  the  solvability  conditions,  one  can  show  that 

(49) 

yn !  yn  i  pn 

where  oja  and  are  the  complex  group  velocities  in  the  x  and  z  directions. 

Since  the  solutions  of  Eqs.  (47)  and  (48)  for  general  initial 
conditions  are  not  available  yet,  we  consider  next  the  special  case  of 
periodic  streamwise  vortices  and  a  single- frequency  disturbance  that  is 
perfectly  tuned  in  the  spanwise  wavenumber.  The  single- frequency 
assumption  corresponds  to  the  case  of  a  disturbance  generated  by  a 
vibrating  ribbon.  The  second  assumption  demands  that  B  =  Bi  and  that 
the  waves  are  modulated  in  the  streamwise  position  only.  Thus,  we 
consider  the  case  in  which  9An/9ti  =  9An/9z1=  0  and  a  =  0.  Then,  Eqs. 

(47)  and  (48)  can  be  rewritten  as 
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(50) 


=  eh  i  A  i  +  eyhi2A2 
=  eh2A2  +  eyh2iAi  (51) 

where 

A  A  A 

hn  =  ^n/9n2»  h  1 2  =  h 1 2 /g 1 2  >  h2i  =  h2i/g22 

It  follows  from  Eqs.  (36)  -  (41)  that  ^i3  =  -  £23  and 

Sin  =  ?2n  for  n  ^  3,  while  it  follows  from  Eqs.  (A. 7)  -  (A. 12)  that 

★  ★  ★  ★ 

the  adjoint  solutions  are  related  by  =  -  c24  and  Cin  =  c 2n  for 
n  +  4.  Hence,  it  follows  from  Eqs.  (A. 2),  (A. 4),  (A. 5),  and  (A. 6)  that 

gi2  =  g22,  hi  =  h2 ,  h32  =  h2i. 

Thus, 

A  A  A  A 

hi  =  h2  dnd  hi2  =  h2i* 


Therefore,  adding  Eqs.  (50)  and  (51)  yields 
(A2  +  Aj)  =  (ehj  +  evhi2)(A2  +  Aj) 
Subtracting  Eq.  (50)  from  Eq.  (51)  yields 
^  (A2  -  Ai)  =  (chi  -  eyhi2)(A2  -  Ai) 
The  solutions  of  Eqs.  (52)  and  (53)  are 


A2  +  Ai  -  2c 


iexp[  J 


(chi  +  cyhi2)dx] 


A2  -  Ai  =  2czexp[ J* (chi  -  evhi2)dx]  (55) 

where  Ci  and  c2  are  arbitrary  constants  that  can  be  determined  from  the 
initial  conditions.  Solving  Eqs.  (54)  and  (55)  gives 

//A  A  ^  A  A 

(eh i  +  evhi2)dx]  +  c2exp[  / (ehi  -  eyhi2)dx] 

Ai  =  Ciexp[J*  (chi  +  eyhi2)dx]  -  c2exp[  J* (chi  -  eyhi2)dx] 

(57) 
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Substituting  for  Ai  and  A2  in  Eq.  (31),  using  Eq.  (35),  substituting 


the  results  into  Eq.  (14),  and  recalling  that  0i  =  $,  we  obtain 


u  =Ci  i 


exp[  J"(i«  +  ehjdx  -  igz  -  iut]  x 
[ciexp(evJ*hl2dx)  -  c2exp(-cvJ* hi2dx) ] 

+  C21  exp[  J(ia  +  dh)dx  +  i 6z  -  iwt]  x 
[c1exp(evJ*h12dx)  +  c2  exp(-evJ*hi2dx)]  +  ... 


(58) 


Equations  (A. 2)  and  (A. 5)  show  that,  in  general,  h12  is  a  complex 
number.  Hence,  one  of  the  terms  multi  pi ing  Ci  and  c2  decays  while  the 
other  grows  exponentially  with  distance.  Thus,  the  growth  rate  a  based 
on  u  for  either  the  wave  with  the  positive  or  negative  8  is 


-  -  a- 


+  e[Real(hj)  + 


1 

Cu 


d£l  1 
3xj 


]  +  el  Real (h12) | 


(59) 


because  Cn  =  C21. 

Equation  (59)  shows  that  the  growth  rate  is  the  sum  of  three 
quantities:  -  o^. ,  the  quasi -parallel  growth  rate;  the  term  proportional 
to  e,  the  effect  of  nonparallelism;  and  the  term  proportional  to  ey,  the 
effect  of  the  streamwise  vortices.  Thus,  in  a  given  physical  situation, 
the  relative  influence  of  the  vortices  and  nonparallelism  depends  on 
the  relative  magnitudes  of  e  and  ev.  For  maximum  amplified  waves, 

_  O 

c  =  0(10  ),  while  for  flows  over  concave  surfaces,  ey  can  be  (0.10), 

depending  on  the  radius  of  curvature.  In  such  situations,  the  effects  of 
the  vortices  dominate  the  effects  of  nonparallelism,  and  the  presence 
of  the  vortices  is  a  very  powerful  instability  mechanism.  Hence,  the 
presence  of  this  mechanism  may  not  be  difficult  to  check  experimentally. 
This  requires  an  experiment  in  which  the  amplitudes  and  spanwise  varia¬ 
tions  of  Gortler  vortices  and  Tollmien-Schl ichting  waves  on  a  concave 
wall  are  measured.  The  spanwise  variations  of  the  Gb'rtler  vortices  and 
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the  generated  Tollmien-Schl ichting  waves  can  be  checked  to  see  whether 
the  wavelengths  of  the  latter  are  twice  those  of  the  vortices.  The 
modification  of  the  mean  flow  due  to  the  presence  of  the  vortices  can  be 

A 

determined  and  the  interaction  coefficient  h12  can  be  calculated  from 
Eqs.  (A. 2)  and  (A. 5)  by  quadrature.  Then,  the  variation  of  u  with  x  and 
y  can  be  calculated  from  Eq.  (58)  and  compared  with  the  experimental 
results. 

It  should  be  noted  that  the  present  analysis  is  valid  only  when 
the  amplitude  of  the  Tol Imien-Schl ichting  waves  Cy  is  small  compared 
with  the  amplitude  of  the  vortices  e  .  As  the  Tollmien-Schl ichting 
waves  grow,  one  needs  to  account  for  their  influence  on  the  vortices. 

In  fact,  they  will  generate  streamwise  vortices  having  an  amplitude 
0(e|)  (Klebanoff,  Tidstrom  &  Sargent,  1962;  Benney  &  Lin,  1960),  which 
may  strengthen  or  weeken  the  primary  vortices,  depending  on  their 
phasings.  This  effect  has  not  been  taken  into  account  in  this  paper. 
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appendix  a 


n2 


n3 


^  +  (u»  -  TT>(Vi!2  +  *  tnie*  )Hy 


’  j  ^ n^n2  +  ^n2^n3  +  ^n3^n 

WQ 

? 

■J  [cn,«;,  +  : 

CO  IK-*; 

/*  2i@ 

=  J  +  ^n^n4  +  ~  +  W*^  +  WV]dy  {A‘3) 

•() 

-h  -?ns.c* 

nn  J  j3Xl  ?n, 


r/M  _  2iav  _  i  da  9;nu  9Un 
L'  0  R  '  3Xi  R  !  dx3  9xj  9xx  ^r\l 


+  v°  d^K2  +  C(u°  -  IT4  SxT2-  -  ihn2  37T  +  ly  (V°;n2)]^ 


3? 


3? 

R~;  3xx  "  R  ^n3  3x7  T  ¥°  3y7 


*  C(u,  -  T2-)  ^  i  t.  &-  ♦  V 


7k;,Jdy 


n3 


(A. 4) 


-2h 


1  2 


'[laUiC2 1  +  V3  |^22  +  ly2-  ^22  -  BiW1^21  -  2igU1C23]c?2 


-2h 


2  1 


R 


+  [iaUxC22  +  ^y  ( V x ^2 2 )  “  3iWi^22  *  2i6VxC23 ]C^3 

+  [iaUlC23  +  (26  -  Bx)HiC23  +  V3  |^x  +  i  |Ml  ^22]^4|dy  (A. 5) 

00  r 

=  j  |[iaUiCu  +  Vi  l^2-1-  +  ly1-  C12  -  61W1S11  +  2i6UiCi3]cJ2 


+  [iaUiCj2  +  |y  (V  1^12)  -  61W1C12  +  2i6Vi^i3]c$3  +  [ i all i c  1 3 

+  (28  -  eJWxCxs  +  Vi  |^-  -  i  |^-  C12]cn]dy  (A. 6) 
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* 

where  the  £nm  are  solutions  of  the  following  adjoint  problems: 


io*nx  +  ^  -  “K2  -  *  -  “2  -  ^n2=  0  (A-7) 

-  DC*t  +  OUo^*2  +  i(U0a  -  -  1  (D2  -  a2  -  B’k*3=  0 


-  +  - 

• 

■  R  <D’  -  «!  -  <S>S.  =0 

(A. 8) 

(A. 9) 

'  °C,  '  i6n 

★ 

V  = 

(A. 10) 

★  ★  ★ 

^n2  ^n3  ^n* 

=  0 

at  y  =  0 

(A. 11) 

★ 

5„m  ■*■  0  as 
nm 

y 

(A. 12) 
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